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ABSTRACT 

Comparison  theorems  are  proved  for  systems  of  equations  of  the  form 

ut  =  d1uxx  + 

Vt  =  Vxx  +  9(x,t,u,v)  . 

Here  u  and  v  are  defined  in  R  x  [0,T]  for  some  positive  time  T,  d^ 
and  <*2  are  positive  constants,  and  f  and  g  are  uniformly  Lipschitz 
continuous  functions  defined  for  (x,t)  e  R  x  [0,T],  (u,v)  e  R  x  R.  It  is 
assumed  that  f  is  a  monotone,  increasing  or  decreasing,  function  of  v 
and  g  is  a  monotone,  increasing  or  decreasing,  function  of  u. 
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SIGNIFICANCE  AND  EXPLANATION 


Comparison  theorem  techniques  have  played  a  central  role  in  the  study  of 
scalar,  nonlinear,  parabolic  differential  equations*  These  techniques  have 
proven  less  successful  in  the  study  of  systems  of  equations  for  several 
reasons.  Usually  a  very  strong  monotonicity  condition  must  be  imposed  on  the 
nonlinear  terms  of  the  equations.  This  severely  restricts  the  applicability 
of  the  comparison  theorems.  Furthermore,  there  are  technical  difficulties 
associated  with  unbounded  domains  for  systems  of  equations  which  are  not 
present  for  scalar  equations.  In  this  report  we  demonstrate  how  these 
difficulties  can  be  overcome  for  certain  systems  of  reaction- diffusion 
equations.  These  systems  have  numerous  applications  including  nerve 

conduction  and  mathematical  ecology. , 

\ 

•  \ 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
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COMPARISON  THEOREMS  FOR  REACTION-DIFFUSION  SYSTEMS 
DEFINED  IN  AN  UNBOUNDED  DOMAIN 


David  Terman 

Section  1:  Introduction 

Comparison  theorem  techniques  have  played  a  central  role  in  the  study  of  scalar, 
nonlinear,  parabolic  differential  equations  (see  [1],  [2]  for  example).  These  techniques 
have  proven  less  successful  in  the  study  of  systems  of  equations  for  several  reasons. 
Usually  a  very  strong  monotonicity  condition  must  be  imposed  on  the  nonlinear  terms  of  the 
equations.  This  severely  restricts  the  applicability  of  the  comparison  theorems* 
Furthermore,  there  are  technical  difficulties  associated  with  unbounded  domains  for  systems 
of  equations  which  are  not  present  for  scalar  equations.  In  this  report  we  demonstrate  how 
these  difficulties  can  be  overcome  for  certain  systems  of  reaction-diffusion  equations. 
Since  the  main  purpose  of  this  report  is  to  demonstrate  how  to  treat  unbounded  domains  we 
do  not  attempt  to  present  the  main  theorem  in  its  most  general  form.  In  fact,  to  better 
motivate  the  results  we  only  consider  systems  of  the  form: 


(1.1) 


Ut  *  d1uxx  *  f(X't'U'V> 

Vt  *  Vxx  *  ***'*'*>  • 


In  a  later  report  we  shall  show  how  to  extend  the  results  presented  here  to  more  general 
systems. 

We  assume  throughout  that  (x,t)  e  where  =  R  *  (0,T)  for  some  positive  T. 

The  functions  f  and  g  are  assumed  to  be  uniformly  Lipschitz  continuous  functions 
defined  for  (x,t)  e  5T  and  (u,v)  €  R  *  R.  The  constants  dj  and  dj  are  assumed  to  be 
positive.  We  must  also  assume  that  f  is  either  an  increasing  or  decreasing  function  of 
v,  and  g  is  either  an  increasing  or  decreasing  function  of  u.  We  shall,  therefore, 
consider  the  following  three  cases: 
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(1.2) 


(a)  f  is  an  increasing  function  of  v  and  g  is  an  increasing  function  of  ul 


(b)  f  is  a  decreasing  function  of  v  and  g  is  a  decreasing  function  of  u, 

(c)  f  is  a  decreasing  function  of  v  and  g  is  an  increasing  function  of  u. 

Systems  of  equations  which  satisfy  one  of  the  above  conditions  arise  quite  often  in 

the  Physical  Sciences.  For  example,  ecological  models  for  two  interacting  species  are 
often  of  this  form.  In  fact,  systems  which  satisfy  (1.2a),  (1.2b),  or  (1.2c)  are  sometimes 
referred  to  as,  respectively,  models  for  symbiosis,  competition,  or  predator  prey.  The 
FitzHugh-Nagumo  equations,  which  model  the  conduction  of  electrical  impulses  in  the  nerve 
axon,  is  another  example  of  a  system  of  equations  which  satisfies  (1.2c).  (See  [4],  (5])« 


The  comparison  theorems  are  presented  in  the  spirit  of  sub  and  supersolutions.  We 
assume  throughout  that  there  exist  functions  (u,v),  (vi#vj#  and  (u,v)  which  satisfy  the 
following  conditional 

(a)  u,  v,  u,  v,  u,  and  v  are  all  bounded  continuous  functions  in  8 , 

(1.3) 

(b)  the  first  derivative  with  respect  to  t  and  second  derivative  with  respect 
to  x  of  u,  v,  u,  v,  u,  and  v  are  all  bounded  continuous  functions  in 

V 

T 

(c)  (u(x,0),v(x,0) )  4  (u( x,0 ) , v( x,0 )  4  (u(x,0 ) , v(x,0 ) )  for  x  e  R, 

(d)  (u(x,t),v(x,t) )  <  ( u( x , t ) , v( x, t ) )  in  8^. 

Here  (a,b)  4  (c,d)  means  that  a  4  c  and  b  4  d.  To  simplify  the  statement  of  the 
comparison  theorems  we  define  the  following  operators: 


(1.4) 


F(x,t,u,v)  5  ut  -  -  f(x,t,u,v) 

G( x, t,u, v)  S  vt  -  -  g(x,t,u, v)  . 


We  are  now  ready  to  state  the  main  theorems. 

Theorem  1 .  Assume  that  (1.2a)  and  (1.3)  are  both  satisfied.  If 

(F(x,t,u,v),  G(x, t,u, v) )  4  (F(x,t,u, v),  G(x, t,u, v) ) 
in  8^,  then  (u,v)  4  (u,v)  in  8^. 
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Note  that  this  theorea  also  implies  that  If  (1.2a)  and  (1*3)  are  satisfied,  and 

(F(x, t,u,v) ,  G(x,t,u,v) )  <  F(x,t,u,v),  G(x,t,u,v)) 

in  Q  ,  then  (u,v)  <  (u,v)  in  Q  . 

T  T 

Theorem  2 >  Assume  that  (1.2b)  and  (1*3)  are  both  satisfied.  If 

F(x,t,u,v)  <  F(x,t,u,v)  and  G(x,t,u,v)  <  G(x,t,u,v) 

in  0  #  then  u  <  u  and  v  <  v  in  0.. 

T  —  T 

Of  course  this  theorea  also  implies  that  if  (1.2b)  and  (1.3)  are  satisfied,  and 

F(x,t,u,v)  <  F(x,t,u,v)  and  G(x,t,u,v)  <  G(x,t,u,v) 

in  Q  ,  then  u  <  u  and  v  <  v  in  Q  . 

T  —  T 

Theorem  3.  Assume  that  (1.2c)  and  (1.3)  are  both  satisfied.  If 

(F(x,t,u,v),  G(x,t,u,v) )  <  (F(x,t,u,v),  G( x, t,u, v) )  <  (F(x,t,u,v),  G(x,t,u,v)) 

in  fl  ,  then  (u,v)  <  (u,v)  <  (u,v)  in  Q  • 

T  — 1  •"  T 

Note  that  in  the  predator-prey  case  (Theorea  3)  one  must  have  a  lower  bound  and  an 
upper  bound  at  the  same  time,  in  the  other  two  cases  one  is  able  to  obtain  one  sided 
estimates.  For  this  reason  the  predator-prey  case  is  usually  the  most  difficult  to  treat. 

A  number  of  authors  have  proved  similar  results  under  the  additional  assumption  that 
the  special  domain  of  u  and  v  is  bounded.  See,  for  example,  Walter  [6] . 

The  stain  tool  in  the  proof  of  Theorea  1  is  the  following  comparison  theorem  for  scalar 
nonlinear  parabolic  equations.  The  proof  of  this  theorea,  with  slight  sKMlifications  since, 
here,  f  and  g  depend  on  (x,t)  aay  be  found  in  [1]. 

Theorea  4.  Suppose  that  u  and  v  are  bounded  continuous  functions  in  Q  and  u. .  v.  . 
.  —  T  t  t 

uxx'  vxx  are  continuous  functions  in  Q^.  Suppose  that  h(x,t,u)  is  a  uniformly 

Lipschitz  continuous  function  defined  in  x  Finally,  suppose  that 
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u  -  u  -  h(x, t,u)  <  v  -  v  -  h(x,t, v)  in  Q_  / 

t  XX  t  XX  T 

u(x,0 )  <  v(x,0 )  in  R  . 

Then  u  <  v  in  ft.* 

T 

Theorems  1 ,  2,  and  3  are  proved  by  constructing  a  sequence  of  functions  (un(x,t), 
vn(x,t)),  n  -  which  are  solutions  of  scalar  equations.  It  is  shown  that  the 

sequence  of  functions  (un,vn)  converge  uniformly  on  bounded  subsets  of  ftT  to  the 
solution  (u,v)«  Hence,  the  proof  of  the  above  theorems  also  gives  us  the  existence  of  a 
solution  to  System  1.1. 

In  this  report  we  only  prove  Theorem  3  since  the  proofs  of  Theorem  1  and  Theorem  2  are 
quite  similar. 


Section  2 :  Proof  of  Theorem  3 

Theorem  3  is  proved  by  approximating  the  solution  of  System  (1.1)  by  a  sequence  of 
functions  (un(x#t),  vn<x,t)),  n  -  1,2,..*,  which  are  solutions  of  scalar  differential 
equations.  We  show,  using  Theorem  4,  that  for  each  n,  (u,v)  <  (u^v^)  *  (u,v)  in  ftp- 
We  then  show  that  the  sequence  of  functions  ^un*vn^  converge  uniformly  on  bounded  subsets 
of  8^  to  the  solution  (u,v). 


The 

functions  un(x,t) 

and  vn(x,t) 

are  defined 

as 

follows.  Let  u^x,!)  be  the 

solution 

of  the  equations 

F (x,t ,u1 ,v)  - 

F(x,t,u,v) 

in 

8 

T 

u1 (x,0 ) 

-  u(x,0 ) 

in 

R 

Let  Vj(x,t)  be  the  solution  of  the  equations: 

G( x, t , u  ,  v  )  *=  G(x,t,u,v)  in  Q  , 

11  T 

v^XjO)  *  v{  x,0 )  in  R 

Assuming  that  the  functions  (Uj ,vt <uk,vk)  have  been  defined  we  let  uk+1(x,t)  be 
the  solution  of  the  equations: 


F(x,t,uk+1 ,vk)  «  F(x,t,u,v)  in  8t  , 
( x,0  )  «*  u( x,0  )  in  R 


We  then  let  vk+i(x*t)  be  the  solution  of  the  equations: 

G(x,t,uk+1,vk+1 )  «  G(x,t,u,v)  in  8t  , 


vk+1 (x,0 )  »  v(x,0) 


in  R 


In  what  follows  it  will  be  convenient  to  set  <u0<x,t),  vQ(x,t))  s  (u(x,t),  v(x,t)) 
and  v^^(x,t))  *  (u(x,t),  v(x,t)).  We  show,  using  induction,  that  for  n  >  1, 


(2.1) 

and 


("1)IVMun-1(x't)  *  (-U^VX't)  «  <-1)n+1un_2(x,t)  in  8t  , 


12*2)  ("1)n+1vn-1(xft)  <  (-1 ) n+1vn(x,t)  <  (-1 )n+1v^_2(x,t)  in  flT  . 
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This 


First  suppose  that  n  *  1.  We  wish  to  show  that  u(x,t)  <  u^x^t)  <  u(x,t)  in  Q^. 

is  proven  using  Theorem  4.  Note  that,  by  assumption,  u(x,0)  <  u^(x,0)  <  u(x,0)  in  R, 

and  F(x,t,u,jO  >  F(x,t,u,v)  *  F(x«t|U^v)  in  Theorem  4  now  implies  that  u  >  u1  in 

ftp*  Moreover,  since  f  is  a  decreasing  function  of  v,  it  follows  that  F(x,t,u,v)  < 

<  F(x,t,u,v)  in  ft^  Hence,  F(x,t,u,v)  <  F(x,t,u,v)  <  F(x,t,u,v)  *  F(x,t,Uj*v).  Theorem 

4  now  implies  that  u  <  u„  in  ft  • 

—  1  T 

Now  suppose  that  (2.1)  holds  for  some  n  >  1.  We  show  that  (2.2)  holds.  If  n  is 
even  then 

u  0(x,t)  <  u  <x,t)  <  u  ,(x,t)  in  ft  . 
n-2  n  n- 1  T 

Since  g(x,t,u,v)  is  an  increasing  function  of  u,  it  follows  that 

G(x, t,u  ,,v  )  4  G(x, t,u  , v  )  <  G ( x, t ,u  ,v  )  in  ft. 
n— 1  n  n  n  n-2  n  T 

Furthermore,  since  G(x,t,u,v)  «  G(x,t,un_j,vn_j)  *  GfXjtjU^v^  *  G(x,t,un_2, vn-2^  in 

Q_.  it  follows  that 
T 

G(x,t,un_2, vn-2)  <  G ( x , t , un_2 , vn ) 

and 

G(x,t,u  , v  )  <  G(x,t,u  .,v  )  in  ft  . 

n- 1  n  n—  l  n—  i  T 

Using  Theorem  4  and  the  assumptions  that  vn_2(x,0)  *  vn-^x»°)  *  vn<x»°)  in  R  we  now 
conclude  that 

v  _(x,t)  <  v  (x,t)  4  v  ( x, t )  in  ft  . 

n-2  n  n-1  T 

Hence,  (2.2)  is  satisfied.  A  similar  arguirtni  shows  that  (2.2)  holds  if  n  is  odd. 

We  now  assume  that  (2.2)  holds  and  show  that  (2.1)  holds  with  n  replaced  by  n  +  1 • 

This  will  complete  the  induction  argument.  If  n  is  even,  then 

v  _ (x, t )  <  v  ( x, t )  <  **  ,(*»*)  in  A-  • 

n—  2  n  n—  i  T 

Since  f(x,t,u,v)  is  a  decreasing  function  of  v  it  follows  that 

r,x't'Vl'vn-2)  <  F,x't'Vl'Vn)  <  F(x't'Un+1'Vll  in  °T  * 

Moreover,  since  F(x,t,u,v)  *  F( x,t,un+1 , vn>  *  F(x,t,un, vn^1 )  -  F(x, tfun_ vn->2)  tn 

P<x't*Vl'Vn-2)  <  Plx't'Vl’V2) 
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it  follows  that 


and 

F(x't'Vi'Vi1  ”  F(x'  'VVi*  in  °T  • 

Theorem  4  and  the  assuaption  that  un<f1(x#0)  «  u^XjO)  *  u^^X/O)  in  K  now  imply  that 

V*'0  *  Un+1(x't)  4  Vl(X't)  in  °T  • 

Hence ,  (2.1)  holds  with  n  replaced  by  n  ♦  1  if  n  is  even.  A  similar  argument  shows 
that  this  is  true  if  n  is  odd. 

We  have  now  shown  that 


and 


"  4  U2  *  U4  *  ***  *  u2n  * 


&  <  v2  <  v4  <  ...  <  v2n  < 


4  U2n+1  *  •••  <  «3  <  »,  <  « 


<  V_  <  •••  <  V_  <  V.  <  V 

2n+1  3  1 


in  flT«  Hence,  there  exist  pairs  of  functions  (U,V)  snd  (U,V)  such  that  t*2ni*2n* 
converges  to  (U,V)  and  ^u2n+1# v2n+1 ^  converges  to  (U,V)  uniformly  on  bounded  subsets 
of  flT  as  n  ♦  •.  Clearly,  (u,v)  <  <U,V)  <  (U,V)  <  (u,v)  in  To  complete  the  proof 

of  Theorem  3  we  show  that  (u,v)  5  (U,V)  5  (U,V)  in 

Let  K(x,t)  be  the  fundamental  solution  for  the  heat  equation.  That  is. 


K(x,t) 


-x2/4t 


OT  * 

-  Vo  -  v> 

Then,  setting  X^  *  d^  '2  and  X^  -  d^  2  ,  we  have  for  each  n  >  1, 
u2n(x,t)  -  /--K(X1x-C,t)u(C,0)dC  ♦ 

♦  J*  /^ltaiX-{,t-T)|f(t,t,BJn,TJ|i.1)  -  F(C,t,u,v)]dCdT 

<2.3) 

v2n(x,t>  *  /^K(X2x-C,t)v(C,0)dC  ♦ 

+  Jj  /^,K(*2x-e,t-T)[g(e,t,u2n,v2n)  -  c(c,t,u,y)]dedt  . 
Passing  to  the  limit,  n  *  m,  in  (2.3)  it  follows  that 
(2.4a)  U(x,t)  -  /^.KJX^C.tluJC.OJdC 

♦  / *  /  jt(X1x-C,t-t)(f(C,t,0,V)  -  F(C,T,u,v)]dCdT 
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(2.4b)  V(x,t)  -  f_J K(^2x-5,t)v(C,0 )dt  +■ 

+  /q  /"„ KtAjX-e.t-lHgfC.T.O.V)  -  G(e,x,u,v)]d£dx  . 

Similarly,  U,  V  satisfy  the  equations 

(2.5a)  U(x,t)  -  /^..KtA^x-C.tjutC.O  )d?  + 

+  /J  -  P(t/T,u,v))  d?dx 

(2.5b)  V(x,t)  “  ! _JtUk2x-Z,t)v(Z,0)dZ  * 

+  /^K(^2x-C,t)[g(5,x,0,V)  -  GU,x,u,v)]d£dx  . 

'Subtracting*  (2«5a)  from  (2.4a)  and  (2.5b)  from  (2.4b)  one  finds  that  in 

(2.6a)  U(x.t)  -  U(x,t)  -  /J  V-e,t-T)  (f  (5,t,U,V)  -  f  ( C, T,U,V) ]  d?dT 

(2.6b)  V(x,t)  -  V(x,t)  -  /g  /"jCU2x-e,t-T)[g(e,T,U,V>  -  g(C, x,U,V) ] d€dx  . 

Let  w(x,t)  *  U(x,t)  -  U(x, t)  and  z(x,t)  ~  V(x,t)  -  V(x,t).  Then  from  (2.6)  and  the 

assumptions  on  f  and  g,  there  exist  bounded  functions  3^(x,t),  t  Y^(x,t),  and 

Y-(x,t)  such  that  in  ft  , 

2  T 

(2.7a)  w(x,t)  -  / g  /^K(X1x-C,t-t)(61(e.X)w(5.X)  +  Y, ( 5,  X )z( 5,x) ]dCdi 

(2.7b)  z(?,x)  -  J*  /“<.K(X2x-e,t-X)[e2(€,X)w(«,x)  ♦  Y2(5,x)z(5(T))dtdT  . 

Let 

B  -  sup  { I B  ( x, t ) |  ♦  | B  ( x , t |  +  |Y  (x,t)|  +  |Y  (x# t ) | > 

(x,t)eftT  1212 

and 

p(t)  -  sup  { | w ( x , t ) |  +  | z ( x , t ) | }  . 
x€R 

Then,  adding  (2.7a)  and  (2.7b)  we  conclude  that 
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p(t)  <  2B  /^ptUdi  for  t  e  (0,T)  . 

From  Gronwall's  inequality  it  follows  that  p(t)  =0  in  (0,T).  Hence, 

(2.8)  (JJ,V)  5  (U,V)  in  fl  . 

It  remains  to  show  that  (u,v)  =  in  8^*  Since  this  follows  from  an  argument 

similar  to  the  one  just  given  we  shall  only  sketch  the  proof.  First  note  that  (2.5)  and 

(2.8)  imply  that  <U,V)  satisfy  the  equations 

U(xft)  -  £jK<X  x-€,t)u(€,0>d*  ♦  J*  /*JC(AlX-4,t-T)]f  (^,t,u,V)  -  F(C,T(u,v)]d£dT 

(2.9) 

V(x,t)  =  /_aBK(A2x-^,t)v(C,0  )dC  +  /q  A2x-C,t-T)  [g(C/T,U,V)  -  G(£,  t  ,u,v)  )d^dx  . 

However,  from  (1.4),  it  follows  that  <u,v)  is  also  a  solution  of  the  equations  (2.9). 
Therefore,  the  proof  of  Theorem  3  will  be  completed  once  we  show  that  the  solution  of  the 
equations  (2.9)  is  unique.  Because  this  follows  from  a  Gronwall  type  argument  similar  to 
the  one  just  given  we  do  not  give  the  details. 
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